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Abstract 

In this paper we discussed the self-adjointness of the Maxwell's 
equations with variable coefheients e and /i. Three different Lagrangian 
are attained. By the Legendre transformation, a multisymplectic Bridge's 
(Hamilton) form is obtained. Based on the multisymplectic structure, 
the multisymplectic conservation law of the system is derived and a 
nine-point Preissman multisymplectic scheme which preserve the mul- 
tisymplectic conservation law is given for the Maxwell's equations in 
an inhomogeneous, isotropic and lossless medium. At last a numerical 
example is illustrated. 



1 Introduction 

Transient electromagnetic field problems are important problems in many 
modern technology applications. M. J. Gotay and J. E. Marsden et al\^ 
discussed the Lagrangian density for electromagnetism on a fixed back- 
ground space-time X with metric and discussed the existence of it's Euler- 
Lagrange equation. In this paper we will discuss this problem from the point 
of inverse variational problem. As a complete depiction of the behavior of 
the electromagnetic field. Maxwell's equations have very interesting struc- 
ture, we will discuss its Lagrangian density and give its Multisymplectic 
Bridge's(Hamilton) form. 

* Supported by the Special Funds for Major State Basic Research Projects , G 1999, 
032800 
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The multisymplectic method is used to phrase the dynamics, which is de- 
scribed by a Euler-Lagrange system, in the terms of the /mzte-dimensional 
space of fields at a given event in space-time. Some Euler-Lagarange system 
can be formahzed in a multisymplectic Hamilton form [21 El E] which mani- 
fest some characteristics of original systems, so we discretize multisymplectic 
Hamilton form directly in order to simulate original system intrinsically. 

Determining whether a given nonlinear partial differential equation(PDE) 
is an Euler-Lagrange equation of some variational problem is fundamental 
work in finite-element method and nonlinear wave theory and etc.. In this 
paper we will use the straightforward result of this problem given by Ather- 
ton and Homsy[5] to discuss self-adjointness of Maxwell's equations in dif- 
ferent representations. We show that the same equations in different forms 
shows different potentialness, consequently it can be expressed as different 
variational problems. So we choose an appropriate form of the Maxwell's 
equations to derive it's Lagrangian and then rewrite it in a multisymplectic 
Hamilton form|2liniE] by using the Legendre transformation. 

This paper falls into 4 parts. In Sec. 2, we obtain the necessary and 
sufficient conditions of self-adjointness for a 1-st order form of the Maxwell's 
equations in an inhomogeneous, isotropic and lossless medium. In Sec. 3, 
we rewrite the equations in a 2-ed order form by introducing two vector 
potential functions. A complex multisymplectic Hamilton form and a mul- 
tisymplectic conservation law are given as well as a ordinary conservation 
law. In Sec. 4, as application of the multisymplectic method, a multisym- 
plectic algorithm is constructed and numerical simulation is shown. 

2 Lagrangian Formalisms for Maxwell's Equations 

First we shall introduce some basic concepts in variational problems. Let 
X C M™", an open subset with smooth boundary dX. M is some function 
space. Given a vector differential equation with a vector of dependent vari- 
ables Ui = X2, • • • , Xm) G M, i = 1, 2---, n, and a vector of differential 
operators N^{u), j = 1,2, ■ • -,71, such that 




(1) 
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We assume the Prechet derivative of N exist, 

The arbitrary function (p now belongs to M. The Frechet derivative derives 
a matrix A'^^ with operator elements N'{i,j) which are the derivatives of the 
operators N''' with respect to uj, so 

K = N'{z,j)l^^,. (3) 

It is called as differential operator matrix of operator matrix A^. Then the 
self- adjoint test functional is given as 



S{^,^)= / i^K^dv= / ^iN'{i,j)cl,jdv. (4) 
Jx Jx 

From Vainberg's theorem the condition of self-adjointness can be given as 

5(V', 0) = 5(0, V) + / divPdv, (5) 
Jx 

if we set the boundary condition to be free, i.e. P vanishes on dX, then 

N'ii,j)=N'ij,ir, (6) 

where N'{j,i)* is denoted as the adjoint operator of N'{j,i). 

Surely a certain variational problem can lead to more than one Euler- 
Lagrange equation and different boundary conditions, responsibly the La- 
grangian and the boundary condition are not unique in the inverse problems 
of finding corresponding functionals for a Euler-Lagrange equation, so the 
requirement for free boundary conditions in Ref.[S] is too strict for some 
applications such as Multi-symplectic method [31 13 which needs to verify 
self-adjointness of a given system to find one of its Lagrangians at first step 
of the method. For example, the wave equation such as 

utt - Uxx = 0, (7) 
according to Ref.|5| its self-adjointness condition is 

^^^|sx = 0, where {t,x) £ X,U2 = Uxx,uu, (8) 

dU2 

so it does not satisfy this condition whatever boundary condition is given, 
but it actually has Lagrangian, therefor we just set self-adjointness condition 
only to be condition @ without free boundary condition. 
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Let us consider a kth order functional on maps u : M"*^'^ — > M^, for 
Maxwell's equation, u = [Hi, H2, H3, Ei, E2, E3]. Here M^'^ is a (1+3)- 
dimensional time-space domain, denoted by coordinates Xj, j = 1, • • • ,4. 

is the target space, denoted by coordinates y*, i = 1, • • • ,6, additionally 
we denote the space TyR^ (g)T*R^'^ by fH. The functional is defined by means 
of a Lagrangian L : it is a function on the space R}'^ x x 9^, i.e. L is a 
smooth function defined on the bundle over M^'^ x with fiber over (x, y) 
equal to 91. We denote this bundle by TM^ «)mi,3xr6 T*R^'^. Hence the 
functional is defined as 

C = [ Lix„u\x),ui{x))oj, (9) 

where uj is some volume 4-form on M}'^, u\, = ^"g^ 1 = 1) 2, • • • , /c. We 
define Legendre transformations as 

BT 

Trr ■.= ^{x,,u\x),nl{x)), (10) 



du] 



it leads to covariant Hamiltonian function 7i on the multisymplectic mani- 
fold r*M6 8)Ri,3x]R6 TRi'^, as 

rL{xj,u\x),^'i)=^'iv^^-L{xj,u\x)y^{x))). (11) 

Now we consider the Maxwell's equations in an inhomogeneous, isotropic, 
and lossless medium 

dH 

H— + VxE = -K 

dE (12) 
-V.H + e- = -J, 

E = [Ei,E2,E3Y is the electric field strength, H = [Hi,H2,H3f is the 
magnetic field strength with superscript T denoting transpose, e is the 
permittivity, // is the permeability, both e and /i are scalar functions of 
space and time variables. J is the external electric-current density, and K is 
the external magnetic-current density. System (|12|) can be arranged in the 
form 

g-z = T, (13) 

where Z = Z{x, y, z, t) is the field vector consisting of the components of 
the electric field strength E and the magnetic field strength H as 

Z=[Hi,H2,Hs,EuE2,Esf (14) 
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simplified into Z = [H,E]'^, similarly JF = J^{x,y, z,t) is the source vector 
composed of the external electric current source J and the external magnetic 
current source K. Q is the operator matrix given by 



Vx -edt 
jjLdt Vx 

where each element of is a 3 — rd order matrix such as 



(15) 



Vx 



/ 

I -m\ eifi)dt = e{f,)- I dt I . (16) 

\~dy as ^ 




T ) , (17) 



Because the operator ^ is a linear operator matrix, i.e. each element of Q 
is a linear operator, its differential operator matrix denoted as Q' equals Q 
itself. So the adjoint operator matrix of Q' is as 

(Vx)T edt + et 
-^dt-^jLt (Vx 

It is obvious that the self-adjoint condition ® is that Q'"^ = Q'* , so the 
equations have to satisfy e = /i, and the two coefficients are independent on 
time variable t. The system in vacuum of course satisfies G'"^ = G'* . The 
systems with such self-adjointness condition have one Lagrangian as 

L=]^<H,y >iH>+]^<E,V xE>-^i<H,Et>-<H,J> ^^^^ 
+ < E,K>, 

where /i = /i(x, y, z) = e(x, y, z). How to calculate Lagrangian density form 
Euler-Lagrange equation is to be introduced in the next section. 
Consider another form of operator matrix Q 

- j^.) . (-) 

with the source function T changed in 

^i = -[^J, -^K]^, (20) 

the dependent variable Z is unchanged. But the differential operator matrix 
of Qi is not a self-adjoint operator matrix, since that 

/^ + i(Vx)T -dt \ 

^'=\ d, i? + i(Vx)Tj^^i' (21) 
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where A and B are two 3-rd order matrix as 



^=|-(f). (-i).l, B 



(22) 

obviously if G'l' = G'*, £ and /i are free of space variables, additionally such 
form excludes the condition of e, /x = 0. 

So we get another Lagrangian of system 1)12^ 

L =— <H,V X H > +— < E,\/ X E > - < H,Et > -- < H,J > 
2e 2// £ 

+ - <E,K > . 

(23) 

This Lagrangian seems equivalent to the Lagrangian (|18|) . but in fact the 
above cases actually dealt with two different systems because of different 
restrictions for the coefficients e and /x in their self-adjointness conditions. 

Both the two Lagrangian are of first order functions, and the Legen- 
dre transformations are unnecessary in constructing their multisymplectic 
Hamilton form, because the conjugate momentums are H and E multiplying 
some coefficients, for example in L2 



8U 1 



dV X H 2e 

dV X E ~ 2fx 
dL2 



E, (24) 
H. 



dEt 

In fact, for linear 1-st order system which is self-adjoint, if its coefficients 
satisfy some conditions, then it can be written in multisymplectic Hamilton 
form directly. 



3 Multisymplectic Hamilton Forms for Maxwell's 
Equations 

Introducing two vector functions U and V satisfying Ut = E and Vt = H, 
so the system can be rewritten in form (2) with different operator matrix G, 
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different dependent variable Z and source function T 



Tlie differential operator matrix Q' equals Q itself. If e and ^ are not depen- 
dent on spatial variables, then it is easy to demonstrate that Q' is self-adjoint, 
so the self-adjointness of the form ((T^ with sketch is verified, and the 
potential functional is given by 



F{v,u)= I f\v,u]g{xz)dxdi 

Jn Jo 



In Jo 

ilfi < V,Vtt >+l< V,Vx Ut > +le < U,Uu > 
n ^ ^ ^ 

^ < [/, V X > +[/ J + VK)dv (26) 

<VuVt> -\ <Vu^ > -\e < Ut, Ut > 
n ^ ^ ^ 

+ ^ < ?7t, V X y > +UJ + VK + dwP)dv, 

where < •, • > represents the inner production of vectors. We know the 
Lagrangian is unique for a given Euler-Lagrange equation under difference 
of a term divP. Here for simplicity we consider the Lagrangian as 



L =^/i < Vt, 14 > +^ < y , V X [/ > 

+ < Ut, Ut> -\<Ut,V xV > -UJ - VK, 
correspondingly the generalized conjugate momentums are 



(27) 



(28) 



We get the covariant Hamiltonian by 

dL 

S = <P,Vt> + <Q,Ut> + < T^^^--, V X V > 

ov X V 

dL 
aV X U 

= < P,H > + <Q,E > -^/i <H,H> -^e <E,E> +UJ + VK, 



(29) 
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here Z = [H, E,V,U, P,Q]'^ . So the Maxwell's equations are transformed 
into the following form 

1. 



-V xU = P - fiH, 
-V xV = Q -eE, 



1. 



1. 



-Pt --V X E 
2 



(30) 



-Qt + ^^xH = J, 



Vt 
Ut 



E. 



The above equations can be organized in the universal form 

MZt + lCV X Z = VzS{Z), 



(31) 



where the rotation action V x Z denotes [V x H,V x E,V x V,V x U,V x 
P, V X Qf^, and /CV x Z consists of three components such that 



/CV X Z = JCiZx + }C2Zy + /C3Z2. 



(32) 



So 



M 



/O 







0/0 

Vo / 






-I 










-/ 



0/ 



( 




2 



V 






2 







where / is the identity element belonging to P?^^ . 
symmetric, and 




R2 







2 


0\ 


2 









































0/ 






(33) 


^A is obviously anti- 


/O -1 






= 1 


») 


, (34) 


\0 


0/ 





18x18 



R^ 



and it is anti-symmetric. Additionally the 



ultimately /Cj belongs to R 
rotation operator Vx = -Ri^ + -^2^ + -^3^- 

The representation (|31|) is simplified expression in vector form, its com- 
plete extension is the multisymplectic Hamiltonian system which was first 



8 



introduced by Bridges and Derks in |4j. It has the multisymplectic conser- 
vative law 

—u + V X K = 0, (35) 
ot 

u and K are the presymplectic forms 

<jj = ^dZAMdZ, n = ^dZAlCdZ, (36) 

where dZ = {dZi,dZ2, ■ ■ •)"'", and for two n-order vector functions p = 
(Pi) ■ ■ '^Pn)'^ and q = {qi, ■ ' qn)'^, notation A is defined as 

n 

dpAdq = ^ dpi A dqi. (37) 

i=l 

The system H31|) also has the local multisymplectic conservation law 
1 . . d . 

-V X [dUAdH - dVAdE] + —{dVAdP + dUAdQ) = 0. (38) 

Now we return to the second case introduced in Sec. 2. We can rewrite 
it in the multisymplectic Hamilton form directly 

-I {)) dt\ E dx\ E )~\ -K/fi 

(39) 

here the covariant Hamiltonian is S = < J, H > —j^ < K, E >. It has 
the multisymplectic conservation law 



V X 



-dHAdH + -dEAdE 
e jj. 



d 

+ —{dEAdH)=0. (40) 



4 A Multisymplectic Scheme for the Maxwell's 
Equations 

Maxwell's equations are simulated by many methods. Its first numerical 
simulation is given by Yee|H] in 1966, which was developed in many more 
efficient algorithms. In this section, we adopt a multisymplectic scheme to 
approximate the system. 
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For simplicity, we consider 1 + 1-dimensional multisymplectic form of 
(|30|) . and suppose that the external magnetic-current vanished, so = 0, 
and the rotation operator becomes 

V X Z = TZiZ^. (41) 

We discretize the equations (|3fl|) by using midpoint scheme in both t 
and X directions. 



At 2Ax 



(42) 



At 2Ax 



-1. — H- 1 1 • I 1 . 

At «+2.?+2' 



1^-1 1 . 



At «+2J+2 

Here the subscript i and j denote respectively x and t directions. The above 
scheme is called the Preissman multisymplectic scheme which preserves the 
multisymplectic conservative law in discrete space, such as 



where 



2Ax At 



(44) 



Eliminating the auxiliary variables V, U, P and Q from the Preissman scheme 
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we get a nine-point multisymplectic integrator as 



7^l- 



4Ax 

Hi+2j+2 + Hi+ij+2 - Hi+2j - Hj+ij 

Hi+ij+2 + Hijjf.2 — Hi^ij — H; 



4At 



7^l 



Hi+2j+2 + 2i?"i+2j+l + Hi+2j — Hij^2 — 2i?"ij+i — Hij 



H+3/2- 



+ ^i+1/2 



4Ax 

Ei+2j+2 + -S'i+lj+2 — -£^i+2j — 



4At 



+3/2j+3/2 + -'i+3/2j+l/2 + -'j+l/2j+3/2 + Ji+l/2j+l/2) f^- 

The numerical process can be depicted by the following expression 

AZ^+^ = BZ^+^ + CZ^ + J, 

where 



A 



C 



D 


Ao 





0... 








o\ 


Bo 


D 


Ai 


0... 














Bi 


D 


A2-- 




















0... 




D 


As 











0... 





Bs 


D 














J 


-D 


Ao 





0... 











Bo 


-D 


Ai 


0... 














Bi 


-D 


A2-- 




















0... 


Bg-i 


-D 


As 











0... 





Bs 


-D 
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/-2D 



B 



\ 



-2D 



J 



( a(0) + 6(0) 
a(l) + 6(l) 



(49) 



-2D) 

where a(i) = J^+i^+i + •^i+fj+i, K*) = + •^*+|j+f ' ^nd 



/ A^7^l 
A^7^l A^7^l 

-AtTei A^7^l 



V 



Set r = 1, 2, • • • , then for every z = 1, 2, ■ 
/eo + ei ei ■ 

ei 61 + 62 62 • 

62 62 + 63 63 



(50) 



-A^7^l / 

, s, • • • , there are 









O--- 6^-2 6r.-2 + 6r._i 6^-1 

V • • • 6r-l 6r_l + 6r/ 

(51) 

here ej = Axe^_,_-^,_(_il G i?^^'^, j = 1, 2, • • • , r. i?, is defined similarly and 
just change ej in Uj = Axnj_^_^^_^il € R^^"^. 

To check the integrator H45() . we consider a simple example. We shall 
take e = /i = 1 and J = 0. Furthermore, we take exact initial conditions as 



(52) 



Hi{x,0)=H3{x,0) = 0, H2ix,0) = -^-sin(x), 
Ei{x,0) = E2{x,0) = 0, E3{x,0) = sin(x), 
under this initial condition we can get a exact solution as following 



Hi — — 0, H2 
Ei = E2 = 0, 



— sin(x — \ / — t), 



£■3 = sin(x - a/^*); 



(53) 
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We compute H2 in rectangles [0, 27r+3] x [0, 0.1], [0, 27r+3] x [0, 1] and choose 
Ax = (27r + 3)/61,At = 0.01. Fig.l and Fig. 2 shows the numerical solutions 
respectively. We can see the solution is moving to x-direction positively. In 
the rectangle [0, 27r + 3] x [0, 10], the trend of the solution is more manifest. 

The absolute error at t-number=1000 is shown in Fig. 4 by the 'error- 
line', and the 'exact solution-line' is the graph of solution 0.01 x H2- The 
former has a very regular route which is so like the solution graph, and their 
vertexes are almost at the zero points of each other. The regular similitude 
give a challenge to modify the scheme which is expected in further task. 

For the multisymplectic form (|39j) there are many existing schemes to 
simulate the system and to preserve the multisymplectic conservation law 
H40|) . such as the usual midpoint scheme which doesn't present here. For 
cases of more general e and fi, the discussion for their multisymplecticity is 
similar, so omitted here. 
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